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Tau-pair process has been analyzed in the ILD detector model as a benchmark process 
for Lol. Results of background rejection, forward-backward asymmetry and polarization 
measurements are obtained with full detector simulation. 

1 Goals for Lol 

Tau-pair process (e+e^ ^Z*,7 r+r^) at ^/s = 500 GeV is one of the benchmark 
processes[T] proposed by Research Director. According to the report, this process is a good 
sample to examine detector performances of 

• tau reconstruction, aspects of particle flow, 

• ttq reconstruction, 

• tracking of very close- by tracks. 

In this process, tau leptons are highly boosted (7 ^ 140), thus decay daughters (mainly 
charged and neutral pions, muons and electrons) are concentrated in a very narrow angle. 
Reconstruction of ttq from two photons is especially challenging for the ILC detectors. 

Observables for the Lol are cross section, forward-backward asymmetry and polarization 
of tau leptons. The polarization measurement requires identification of tau decays, including 
reconstruction of ttq. Efficiency and purity of event selection cuts should be also a good 
measure of detector performance. 

For physics motivation, tau-pair process is important as a precision measurement of the 
electroweak theory. For example, measuring cross section and forward-backward asymmetry 
of tau-pair process very precisely can probe existence of heavy Z' boson. 

2 Analysis framework and events 
2.1 Event samples 

Events of ILD_00 Lol mass production[2] are used for this study. Events reconstructed and 
listed at DESY by approximately end of February are used in this analysis. 10.3 M SM 
events generated in SLAG are processed for background estimation with appropriate event 
weight. 

Since the SLAG events have a polarization issues for tau-pair events, tau-pair events 
generated in DESY are used instead of SLAG events in this analysis. For other modes 
including tau, SLAG events are used. Whizard 1.51 and TAU0LA[3j are used to generate 
the DESY events. Statistics of the signal channel is 500 fb^^ both for e^ej^ and e^e^ (total 
2.3 M events). 
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Bhabha process (e+e~ elastic scattering) is an important background for tau-pair anal- 
ysis. Since the cross section of Bhabha process is too large (~ 17 nb for each polarization 
in SLAC events), following preselection is applied to the SLAC events before simulation. 

• I cos6'| of electron or positron must be smaller than 0.96. 

• Opening angle between electron and positron must be larger than 165 deg. 

After the preselection, the cross section is reduced to 50-90 pb. ~1.0 fb^^ of preselected 
Bhabha events are simulated. 

Preselection is also applied to 77 — > tt events with following cuts: 

• Opening angle between two taus must be larger than 170 deg. 

• Energy sum of two taus is greater than 30 GeV. 

The total cross section after the cuts is around 18 pb. Around 100 k events passing prese- 
lection are processed. 

Integrated luminosity is assumed to be 500 fb~^ each for two polarization setups, c^c^ 
and ej^e^. Assumed polarization ratio is 80% for electron and 30% for positron (i.e. for 
e£e^ setup 90% of electrons are leftly polarized and 65% of positrons are rightly polarized). 

2.2 Tau clustering 

For tau clustering, an original clustering processor (TaJet) is applied to the output of Pan- 
doraPFA. Following is a procedure of the processor. 

1. Sort particles in energy order. 

2. Select the most energetic charged particle (a tau candidate). 

3. Search particles to be associated to the tau candidate. Criteria is: 

(a) Opening angle to the tau candidate is smaller than 50 mrad., or 

(b) Opening angle to the tau candidate is not larger than 1 rad. and invariant mass 
with the tau candidate is less than 2 GeV (m^ = 1.777 GeV). 

4. Combine energy and momentum of the tau candidate and associated particle and treat 
the combined particle as the new tau candidate. 

5. Repeat from 3. 

6. After all remaining particles do not meet the criteria, remaining most energetic charged 
particle is the next tau candidate. (Repeat from 2.) 

7. After all charged particles are associated to tau candidates, remaining neutral particles 
are independently included in the cluster list as neutral fragments. 

In the clustering stage, events with > 6 tracks are pre-cut to accelerate clustering since 
> 99% of tau decays have < 3 charged particles. Event with only one positive and one 
negative tau clusters are processed with latter analysis. 
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Table 1: Cut statistics for background suppression. Preselection (See Section [^TT] for details) 
is applied for Bhabha events before these cuts. Number of events are normalized to 500 
fb^^. The same statistics is used for (a) and (b): only event weighting is different. 



3 Background suppression 

Main background of tau-pair analysis is Bhabha (c+e~ — >e+e~), WW Ivlv and 77 — >• 
T+r". Since cross sections of Bhabha and two-photon events are huge (about 10^ and 10'^ 
larger than signal, respectively), we need tight selection cuts for those background events. 
Following cuts are applied to signals and all SM background events after the tau clustering. 

1. Number of tracks < 6. Included as a pre-cut in tau clustering processor. 

2. Only one positive and one negative tau clusters must exist in the event. (Neutral 
clusters are allowed.) 

3. Opening angle of two tau candidates must be > 178 deg. 
This cut efficiently suppresses WW — » Iviv background. 

4. ee and /i^ events are rejected. 

Charged particles depositing > 90% of their energy in ECAL are identified as electrons, 
and charged particles depositing < 70% of their energy (estimated by curvature of their 
tracks) in ECAL-I-HCAL are identified as muons. Events with two electrons or two 
muons are rejected in this cut. This cut is especially for suppressing Bhabha and 
e+e^ l^.'^IJ.^ events. Signal loss is about 6%. 

5. I cos6'| < 0.95 for both tau clusters. 

t-channcl Bhabha events arc almost completely suppressed by this cut. 20% of signal 
events are lost. 

6. 70 < Evis < 450 GeV. Evis does not include energy of neutral clusters. 

Lower bound suppresses 77 t^t^ events, and upper bound suppresses Bhabha 
events. Signal lost is negligibly small. 
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Figure 1: Distribution of cut values. Left column shows c^g^ distribution and Right column 
shows e^e^ distribution. Cuts are applied from top, and ec and jiji veto cuts are applied 
between second and third rows, whose distributions are omitted. 
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Figure 2: Angular distribution of t+ momentum direction. Number of events are normalized 
to 500 fb^^. Error bars stand for statistical errors in current MC statistics. The same 
statistics is used for (a) and (b): only event weighting is different. 



Table[T]shows the result of these cuts and Figure [T] shows distribution of cut values. Most 
of the background is effectively cut off by the cuts, ^^10% level of the signal. Remaining 
background is mainly Bhabha, 77 tt and WW^ £i^£iy. 

Purity of tau selection is 92.4% in e^e^ sample and 93.6% in epe^ ■ The difference is 
mainly from difference of the cross section between each polarization. 

Selection efficiency of tau-pair events is literally low (15.8% in e^e^ and 16.3% in Cj^e^ )• 
However, the 'nocut' number contains radiative events, which have effectively lower ^/s and 
should not be used in the analysis. These radiative events are cut by the opening angle 
selection. The real efficiency varies by the definition of the events. The acceptance of 
'softly-radiated' tau-pair events is determined by the opening angle cuts. Loosing the cut 
accepts more events, although 77 tt and WW^ £i^£iy background significantly increase. 



4 Cross section 

Cross section can be easily obtained by count-based method since background amount is 
low. Assuming background subtraction can be performed in the error of statistics, we obtain 
number of signal event as 125400±368 (eZe^) and 103197±332 (cj^e^), ie. 0.29% and 0.32% 
statistical error, respectively. The statistical error is dominated by signal statistics, so poor 
statistics of background events in the current MC sample can only have small effect on these 
numbers (0.30% and 0.33% statistical error, even if background is doubled). Systematic 
error can be introduced by polarization error, MC incorrespondance to real detector etc., 
but it cannot be accurately estimated in this stage of detector development and thus not 
considered now. 



5 Forward-backward asymmetry 

Figure [2] shows a result on angular distribution of t+ Icptons (t^ events are essentially the 
same event- by- event since we require opening angle > 178 deg.). 
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Assuming that background can be subtracted effectively, forward-backward asymmetry 
is calculated by following formulae. 



-^^B = l/(^-iV.) (2) 




aNp = VNf + Nfbg, (^Nb = VNb + Nbbg, (3) 



where Nb is number of signal events in backward region (cos0 < 0), Np is number of signal 
events in forward region (cos0 > 0), Nbbg is number of background events in the backward 
region and NpBG is number of background events in the forward region. The formulae can 
be reduced to 

2^Nl{Np + Npbg) + N],iNB + Nbbg) 

""^^^ ^ {Np + Nb)^ ■ 

Result of the calculation is: 

eLe,^ : A^f = 95529, iVs = 29872, iVi^sG = 9201, TVbbg = 1130, Afb = 52.36 ± 0.259^5) 
e^e^ : A^f = 75556, iVs = 27640, TV^sG = 5477, iVsBG = 1605, A^s = 44.19 ± 0.289^6) 

Statistical accuracy of ApB is 0.48% and 0.63%, respectively. 



6 Decay mode separation 

Separating decay modes of tau is essential for the polarization measurement. There are five 
dominant decay modes of tau, t+ e^v^Vr (17.9%), t+ ^'^TT^Vr (17.4%), t+ tt^Vt 

(10.9%), T+ p+Vr TT+Tl^Vr (25.2%), and T+ -A+Vr ^ TTTTTTZ^^ (9.3% (l-prOUg) 

and 9.0% (3-prong)). Other decay modes (10.3%) include Kaons and multi-pions in other 
resonant modes or continuum. 

We utilize a neural network for the decay mode selection. Two separate networks are 
trained for 1-prong and 3-prong events. 1-prong decay includes e^V^Vr, fJ-'^i^i^r, '"""''^'r, P^i^t 
and Sk^i/r modes. Input variables of the 1-prong neural net are as follows. 

• Two lepton-ID values. Likelihood-based lepton ID software was developed, but due 
to the known issues of the event production the lepton ID is not properly worked on 
the mass production sample. As a simpler lepton ID, we use ratio between the energy 
deposit of the electromagnetic calorimeter (ECAL) and the total deposit energy for the 
electron ID, and ratio between the calorimeter energy deposit and the track momentum 
for the muon ID. These two variables are included in the neural network. (Variable 
(a) and (b).) 

• Energy of the charged particle and two kinds of energy sums of the neutral particles. 
The neutral energy sums contain particles whose ECAL energy deposit is > 80%, and 

< 80% of the total energy deposit, respectively. Particles with ECAL energy deposit 

< 80% are considered to be hadrons, which contain more spurious particles from a 
mis- fragmentation of energetic charged particles mainly at HCAL. The energy sums 
are especially expected to discriminate tt mode. (Variable (c), (d) and (e).) 
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• Number of neutral particles except neutral hadrons. Number of photons is a powerful 
information to separate p (expected number of photons is 2) and ai (expected number 
of photons is 4). (Variable (f).) 

• Invariant masses of all reconstructed particles except neutral hadrons and invariant 
masses of photons. Invariant masses of all reconstructed particles should equal to 
the masses of intermediate particles, p and ai. If photons are reconstructed properly, 
invariant masses of photons are close to that of ttq. For the photon / hadron separation, 
above criteria is used again. (Variable (g) and (h).) 

• Energy of the third-energetic neutral particle. This variable is also to separate p and 
ai. Since p can have at most two photons, energy of the third photon should be small 
even if it exists in the p mode. (Variable (i).) 

We use two hidden layers, first layer has 18 neurons and second has 9 neurons. The 
output neurons are likelihood value of e^v^Vr, P-^vJI^t: tt^j^t, P^Vt and Sk^Vr modes (5 
neurons), which is set to 1 (true) /O (false) by the MC information in the training samples. 

For the 3-prong events, only ai is the discriminating decay mode. Input variables (a), 
(b), (c), (d), (f), (g), (h) (noted with (a')-(h') in Fig. 3]) in the 1-prong selection are also 
included in the 3-prong selection. There is one additional variable, which is: 

• Invariant mass of all charged particles. This should equals to the mass of ai if the 
decay is ai mode. (Variable (j').) 

We use two hidden layers with 16 and 5 neurons. The only output neuron stands for 
likelihood value of ai mode, set to 1/0 in the training samples as well. 

For the training, half of the tau-pair events in the mass production are used. Number of 
epochs is 1000 for both 1-prong and 3-prong network. 

Fig. [3] and S] shows distributions of the input variables, and Fig. [S] shows distribution 
of the output neurons. The mode selection is applied based on the values of the output 
neurons, as follows. 

• If one or more of the values of the output neurons exceed 0.5, The neurons which gives 
the highest output value is used as the selection. 

• If no output values exceed 0.5, the event is classified as 'others'. 

Table [5] shows the obtained efficiency and purity for the mode selection. These values 
are obtained with the half of the tau-pair events which arc not used for the training. > 90% 
efficiency and purity is obtained for all decay modes except 1-prong ai decay. 

7 Polarization Measurement 
7.1 Optimal Observable 

To identify tau polarization, optimal observables[3] are used for e^V^Vri P^i^i^r, tt+i^t- and 
p+i^i- decay modes. Decay distribution of all tau decay can be described as the same form, 

W = ]^{l+pcoseh) (7) 
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Figure 3: Distributions of the input variables for the 1-prong neural network, e^ (80%) 
polarization is used for the plots. 





Figure 4: Distributions of the input variables for the 3-prong neural network, 
polarization is used for the plots. 
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Figure 5: Output variables for the neural net selection, (a)-(e) are the output of the 1-prong 
neural net for e'^V^i^r, fJ-^i^i^T, t^^^t, P^^t and a^Vr rnodes, respectively, (f) is the output 
of the 3-prong neural net for &^Vt identification. Cj^ (80%) Cp polarization is used for the 
plots. 
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Table 2: Purity and efficiency of the tau decay mode selection with neural networks. Process 
background is not included in the purity & efficiency numbers. 
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where p is polarization of r (-1 to 1) and 9^ is the opening angle of polarimator vector h with 
respect to the r momentum vector. Explicit notation of h varies by the decay modes: for 
pure-leptonic decay, flight direction of antineutrino is h and for ■k'^Vt decay, flight direction 
of pion is h. For the multipion decay, h is constructed from the hadronic current. 

To reconstruct p from a set of observables, we split W to p-dependent and p-independent 
components such as 

W{0 = f{0+P9{0. (8) 
and the optimal observable lj is defined as 

.^m. (9) 

/(e) 

By definition, probability density P at w for the polarization p gives 

P(w;p)-P(w;p = 0) 



P{uj;p = 0) 



uj (10) 



and p can be easily obtained from the lo distribution. 

The explicit formula of lo for each decay mode is as follows [5]. 

1. Pure-leptonic decay: 

where x is the lepton energy divided by r energy (250 GeV in this case). Since the 
pure-leptonic decay mode has two missing neutrinos, polarization discrimination power 
is weaker than semi-leptonic decay modes. 

2. TT^Vr decay: 

w^ = 2x-l. (12) 

This mode has maximum polarization discrimination power since h can be fully recon- 
structed. 

3. /O^i^r decay: This decay mode has multiple observable particles and thus more compli- 
cated formula to describe uj. Tau momentum direction is unobservable in this decay, 
so it is integrated out in the tu formulation. The explicit formula is: 



V'-l 3coB-^ 0-1 \ . , / 3coB^ 0-1 , 



(13) 

cosih = — (14) 

{ml-Q^)^x^ -AQ^/s 

X = 2^ (15) 

where is the energy sum of p (which equals to the cluster energy), is the 
invariant mass of the visible particles (should equals to rup = 0.77 GeV but obtained 
from the event), y/s is the center-of-mass energy (500 GeV), is the angle of the p 
flight direction with repect to r direction in r-rest frame, and (3 is the angle of the 
charged pion flight direction with repect to p direction in p-rest frame. 
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7.2 Polarization measurement 



Figure [6] shows the uj distribution for each decay mode passing the neural net selection. For 
the Icptonic mode, most of the events are concentrated on the uj '-^ Q region, reflecting to 
the weak discrimination power. For the 'k'^Vt and p^v^ modes, distribution is broadly 
distributed and large difference between left and right polarization can be seen. 
Polarization p can be obtained by following procedure. 

1. Mode and process background is eliminated from each bin of the w histograms and 
statistical error of background remains included in the error of each bin. 

2. Histograms from all decay modes are summed into one histogram. 

3. The histogram with polarizing sample (polarization p) is divided by non-polarizing 
sample after normalizing both histograms. 

4. Perform linear fit passing (0,1) to the divided histogram (one parameter fit). Obtained 
slope stands for -p. 

Figure [7] shows the combined lo distribution and Figure [5] shows the linear fits to obtain 
V value. Obtained p is -63.82±0.66% (eLe,^, 80% and 30%) and 50.83±0.79% (e]^e+, 80% 
and 30%). 

8 Summary 

Tau-pair process has been analysed in the ILD_00 detector model. After the tau selection 
cuts, statistical error of cross section measurement is 0.29% (clCj!;, with 80% and 30% 
polarization, respectively) and 0.32% (cj^e^). Process background can be suppressed to 
around 10% of signal events. Forward-backward asymmetry can be determined with 0.48% 
and 0.63% statistical error. 

Polarization measurement needs separation of decay modes. The neural net selection 
gives > 91% efficiency and > 88% purity of mode selection for all major decay modes except 
a\v \ prong mode. Polarization analysis of c+I^^'t-, fi^T^i'r, ti^Vt and p^Vr decay mode is 
performed using the optimal observable method, and it results in P{t) = —63.82 ± 0.66% 
ie£e+) and P{t) = 50.83 ± 0.79% ieZe+). 

The SLii>r mode is not included in the current polarization measurement. For the 3-prong 
ai decay, r direction can be reconstructed from the vertex information and it can improve 
the analysis power to the same level as tt^Vt mode. However, since the branching ratio 
of 3-prong ai decay is only about 9%, the expected improvement with 3-prong ai decay is 
about 20%. 
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Figure 7: Distribution of the optimal observable after summing up all decay modes. 
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Figure 8: Ratio of the polarizing sample to the non-polarizing sample. Linear fit is applied 
to obtain p value. 
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